Work and its fluctuations in a driven quantum system 
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We analyze work done on a quantum system driven by a control field. The average work depends 
on the whole dynamics of the system, and is obtained as the integral of the average power operator. 
As a specific example we focus on a superconducting Cooper-pair box forming a two-level system. 
We obtain expressions for the average work and work distribution in a closed system, and discuss 
control field and environment contributions to the average work for an open system. 



The fluctuation relations (FRs) [TJ [5] govern work 
and dissipation in small classical systems when they are 
driven out of equilibrium. They have recently attracted 
lots of attention because of their applications in molec- 
ular systems [3]. Fluctuation relations can also be accu- 
rately studied in single-electron transport jlHB]. A nat- 
ural question is if similar concepts and experiments can 
be extended to the quantum regime. The first attempts 
in this direction focused on finding a proper work opera- 
tor [T|r7Ul0j. However, after a long debate, it has become 
clear that this approach has serious drawbacks [IT] . Work 
is characterized by a process, not only by the instanta- 
neous state of the system [TTJ [T3] , and therefore it can- 
not be defined by an operator local in time, which would 
disregard the actual evolution of the system under the 
driving protocol. Although this is not an issue for closed 
systems it can become critical when discussing work in 
open systems. Alternatively the work has been defined 
through a two- measurement process (TMP) [TUHTB] . The 
energy of the system is measured at the beginning and at 
the end of the evolution and the work done in a process is 
determined by the corresponding energy difference. This 
definition has the advantage that the quantum FRs can 
be immediately obtained and they resemble the classical 
ones. In this proposal the system does not interact with 
the environment and, thus, the dynamics is unitary. 

To circumvent the problem of extending the TMP ap- 
proach to an open system |17l I18j . we introduce work 
in analogy to that in the classical case as an integral of 
the injected power during the evolution. Let the evo- 
lution of the system be governed by a time-dependent 
Hamiltonian H(t) driven by a control parameter X(t). 
The corresponding power operator is then given by P — 
dH/dXX = dH/dt. If the state of the system is described 
by its reduced density operator p(t), the average power 
is given by (P(t)) = Tr{p(t)P(t)} and the expectation 
value of the work done on the quantum system is 

(W) = f T (P(t))dt. (1) 
Jo 

This way, the work explicitly depends on the whole evo- 
lution of the system through p(t) containing the informa- 
tion about the dynamics which can be unitary or not. To 



address this point, we differentiate the average energy of 
the system, (H) = Tr{p(t)H(t)}, yielding 

i{E) = T, { fH } + { P) (2) 

Under quite general assumptions the dynamics of the re- 
duced density operator of the system can be described by 
a master equation [19 dp/dt = —i[H,p] + C(p) where 
the contributions on the right-hand-side are given by the 
unitary and dissipative dynamics, respectively. By sub- 
stituting the above result into Eq. ([2]), we find that 
there is no contribution due to the unitary dynamics 
since Tr{[H,p]H} vanishes. Then the average power 
reads (P(t)) = d (H(t)) jdt - Tr{£(p)H(t)} and the cor- 
responding average work is given by 

(W) = (H(T)) - (H(0)) - [ T dtTr{C(p)H(t)}. (3) 

Jo 

If the system does not interact with the environment, 
only the first difference on the right-hand-side survives 
in Eq. (J3j) , and the average work is equal to the variation 
of the internal energy. The last term describes the en- 
ergy exchange with the environment during the evolution 
process, and is dependant on the particular realization of 
the evolution trajectory. We call this contribution aver- 
age heat and denote it as Q. In thermodynamical terms, 
Eq. ([3]) is the first law in the quantum regime, and it has 
been discussed previously in Refs. [301 HI] as the energy 
balance equation. 

The average work definition in Eq. ^ is more gen- 
eral than the TMP since it takes into account the full 
quantum evolution. For a closed quantum system, given 
the initial quantum state and a driving protocol, the 
evolution is determined completely by the Schrodinger 
equation, while the initial measurement performed in the 
TMP causes a collapse of the quantum state before the 
beginning of the protocol. Thus, the two approaches 
yield different results for the initial states with coherent 
superpositions of the eigenstates of H(0) and the TMP 
result is recovered if the system is initially in an eigen- 
state or in an incoherent superposition of the eigenstates 
of.ff(0). 
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The second advantage of Eq. (|3| is that it allows to 
define and calculate the heat only in terms of the system 
density matrix (and its dynamics) and, thus, it avoids 
the formidable task of measuring the environment de- 
grees of freedom as in the TMP [T21 [13]. In addition, 
no assumption has been made on the dissipative dynam- 
ics then Eq. ([3| can be used with any master equation 
(Lindblad, non-Lindblad, non-Markovian) . 

Cooper-pair box as a driven quantum two-level system. 
We consider a Cooper-pair box (CPB) [2"2"H2"4"| consisting 
of a superconducting island connected to a superconduct- 
ing lead by a Josephson tunnel junction. The system is 
described by the circuit scheme in the inset in Fig. [2] 
and it is characterized by a voltage source V g , coupling 
gate capacitance C g , a Josephson junction with energy 
Ej and capacitance Cj. We denote Cs = C g + Cj. Re- 
sistor R, to be discussed in the last part of the paper, 
forms the dissipative environment of the box. This sys- 
tem is a good candidate for a calorimetric measurement 
of quantum work distribution |25j . 

In the regime e = Ej/(2E C ) < 1, where E c = 2e 2 /C s 
is the charging energy of the box, we can treat the CPB 
as a two-level quantum system. Denoting with |0) and 
|1) the state with zero and one excess Cooper-pairs on 
the island, respectively, the Hamiltonian reads 

H = -E cq (\l)(l\ - |0)(0|) - ^(|1)(0| + |0)(1|), (4) 

where q = C g V g /(2e) — l/2 is the normalized gate voltage. 
We assume driven evolution: a linear gate ramp q(t) — 
— 1/2 + t/T over a period T starting from t = 0. The 
ground and the excited states of the system are separated 
by the energy gap hujQ = 2Ec Vl 2 + e 2 which reaches the 
minimum Ej at t = T/2, see Fig. fll We recover the stan- 
dard Landau-Zener (LZ) model [261127] where the system 
is excited when driven in a non-adiabatic way through a 
avoided crossing. The time-dependent eigenstates of the 
Hamiltonian (|4) are \g) = ~^{V^ 



7?|0) + yTT^|l))and 
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V2 



(y^TTvlO) - \/l~ r ^|l)) where n = q/yjq 2 + ■ 



For this system, the power operator is P = Ecq{^ — 
2n), where n = |1)(1| is the operator of the number of 
Cooper pairs on the island and 1 = |1)(1| + |0)(0| is 
the identity operator. We calculate the time-dependent 
average of the power operator in the Heisenberg pic- 
ture, P H {t) = W(t)PU(t), with the time evolution op- 
erator U(t), and the state \ip(0)) that does not change 
in time. Here we focus on the first and second mo- 
ments of the work done on the CPB, which can be 
expressed through P H {t) as (W) = dt(P H (t)) and 

(W 2 ) = jjdt^jdt^it^it,)), where (...) = 
(-0(0)1... |^(0)). Explicitly: 



(W) = E c (l-l f{n H {t))dt 



(5) 



J (7 Hi) 




FIG. 1: (Color online) Schematic presentation of avoided 
crossing with the eigenstates of energies E g and E e in a CPB 
as a function of the normalized gate charge q. The phases 
which contribute to the interference are explicitly indicated. 



and 



(W 2 



-E 2 C 



= 2E C (W) 
X -T 2 l dh 



(6) 



T 



dh(h H \t 2 )n H \ti)) 



Equations ^ and J6 1 can be applied for both closed and 
open systems [T^J EH . 

Instantaneous transition regime, unitary evolution. If 
the time of the control ramp is much shorter than the re- 
laxation and dephasing times, the evolution of the system 
can be considered unitary. For e <C 1, the LZ transitions 
are localized near the minimum energy gap at t = T/2 
and the dynamics is well approximated by the instan- 
taneous transition model [311 [32], i.e., the evolution is 
composed of pure adiabatic evolution and instantaneous 
LZ transitions at t = T/2, see Fig. [I] All work, spent 
exactly in these LZ transitions, is stored in the system 
(CPB) as increased internal energy. Along the adiabatic 
region, the evolution operator reads Uk{t) = exp [— i£|<7 z ] 
where £| = f t k i dTU)( S {r)/2 is half of the integrated en- 
ergy gap, tk and ifc_i denote the time limits of the adi- 
abatic evolution and cr z = \e(t))(e(t)\ — \g(t)) (g(t)\. The 
transfer matrix for the instantaneous LZ transitions in 
the basis {\g(T /2)), \e{T /2))} reads 



N, 



LZ 



LZ 



-yp , 

Vi-Pi 



LZ 

lz&~ 



(7) 



where <p = 8 {log 5 — 1) + argT(l — iS) — tt/A (T is 
the gamma function) [311 [32]. The probability of the 
LZ transition is given by P^z — e~ 27r<5 , where 8 — 
E 2 T /{8E C ) = e 2 E c T /2. We consider the system ini- 
tially in a superposition of ground and excited state and 
\ipo) = ot\g(0)) + Vl - a 2 exp («7)|e(0)). We can then 
write for t < T/2, U{t) = Ui(t), while after the LZ 
transition at t > T /2, U(t) = U 2 {t)N LZ U 1 {T /2). The 
calculation of \tp(t)) is in the Supplementary Material 
(SM)[33L Here we discuss only the relevant results. 
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FIG. 2: (Color online) (a) Average work (W) normalized to 
Ec for different initial states: a = l/v2 and 7 = (blue os- 
cillating curve), ground state a = 1 and 7 = (top black ex- 
ponentially decaying curve) and excited state a — and 7 = 
(lower black exponentially decaying curve). The purple curve 
denotes the behaviour \J Plz(X — Plz)- (b) Work variance 
{SW 2 } normalized to E c . The blue solid lines contain the os- 
cillatory behaviour for initial state a = l/v2 and 7 = 0. The 
black dashed line is obtained for initial state, a = 1 and 7 = 
(ground state). We have used e = Ej/(2Ec) = 0.05. Inset: 
Circuit scheme of the Cooper-pair box (CPB) connected to a 
dissipative environment R. 



From Eq. ([5]), the corresponding average work is 

(W) = E c l(2a 2 - l)P LZ (8) 

+2ay/l-a 2 y/(l-P LZ )PLZ cos( 7 + <p + 2ff )]. 

The first term represents the work done on the system 
which is initialized in the ground or in the excited state, 
i.e., a = 1 or 0, respectively. The second term with its 
characteristic oscillatory behavior is due to the quantum 
interference at the LZ avoided crossing [32] . This addi- 
tional contribution is always present when the system is 
initially in a coherent state and it is a clear difference 
with the respect to the TMP. This difference is high- 
lighted in Fig. [2] (a) where we plot the analytical result 
111 Eq. (9) as a function of T for different initial states. 
The oscillating behavior of (W) is obtained for a = \j\/2 
and 7 = and we should compare it with the prediction 
of the TMP (W) = 0. The two exponential decays with 
Plz are obtained for ground and excited initial state, 
a = 1 and a — with 7 = 0, respectively. The behav- 
ior for thermalized initial density matrix can be obtained 
from this two curves with the correct weighted average. 
With the same approach, the evaluation of Eq. ^ 



yields (W 2 ) = PlzE c for the second moment indepen- 
dent of the initial state. With these results the corre- 
sponding rms fluctuation of work can be immediately 
calculated as (5W 2 ) = (W 2 ) -(W) 2 . Figure |j(b) shows 
the behavior of (8W 2 ) for different initial states. Numer- 
ical simulations confirm the analytical results presented 
in Figs. [2] 

The definition of work in Eq. ^ and the TMP give 
the same results if the system is initially in an eigenstate 
of H (0) or an incoherent superposition of them. In the 
interesting case in which the system is initialized in the 
ground state, i.e., a — 1, and for nearly adiabatic drive 
{Plz € 1), we have a linear response result linking the 
average work and its fluctuations as (SW 2 ) — Ec(W). 

In this specific case, the first two moments of work 
calculated above agree with the full work distribution 
p(W) which for a closed system with unitary evolution 
U(T) can be found essentially by direct comparison of the 
initial, -ff(O), and final, H(T), Hamiltonian of the sys- 
tem. Indeed, the work generating function G(u) (Fourier 
transform of the distribution) can be written as (see, 
e.g., [Ill): G(u) = T^{U^(T)e iu6 ^U{T)e- iu6 ^p } 
where po is the initial density matrix of the system as- 
sumed to be diagonal together with the initial Hamilto- 
nian -ff(O). For the CPB considered above and the sys- 
tem initialized in the ground state, this equation gives 
G(u) = 1 + PLz{e mEc — 1) which corresponds to the 
following work distribution: 

p(W) = (1 - P L z)S(W) + P LZ 5(W - E c ) ■ (9) 

This distribution agrees with the first two moments and 
can be used to find the higher moments. 

Open system with slow and fast relaxation. The most 
interesting and non-trivial test of Eq. ^ is when the 
system interacts with the environment during the evolu- 
tion of q. To evaluate the heat contribution in Eq. ^ 
we need to consider a concrete example of the system- 
environment interaction. If the time and the ramp time 
are of the same order, dissipation Q takes place dur- 
ing the driven evolution. To evaluate dissipation during 
the sweep, we then solve the master equation (ME) of 
the CPB adapted from the corresponding ME of Refs. 
[251 [51] . This ME and some details of the analysis are 
given in the SM. The environment is described by the 
resistor R coupled capacitively to the island of the CPB 
(inset in Fig. [2]). As above, we assume that the temper- 
ature is low as compared to the excitation energy. If the 
system is initially in the ground state, the average heat 
released to the environment during the ramp normalized 
by the total work done for a few values of e is shown in 
Fig. [3] (a) based on the numerical solution of the ME 
(solid lines). Dependence of the same quantity on the 
minimum energy gap Ej (for scaling purposes the hori- 
zontal axis is e 2 ) is shown in (b). The apparent depen- 
dences on the various parameters in Fig. [3] are captured 
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The system tends to follow the instantaneous equilibrium 
state but due to the drive there are corrections to this dy- 
namics which cause the heat emission. We calculate the 
heat in this regime directly from the semi-classical mas- 
ter equation. As shown in SM, the same result is still 
obtained if one takes into account quantum corrections 
from the finite drive to adiabatic dynamics. 

The integrand in the heat contribution in Eq. ^ can 
be conveniently written as —hp gg ujQ. If we denote with 
P ge and T eg the excitation and relaxation rates, respec- 
tively, in the semi-classical limit we have p gg = —T^p gg + 
T ge where T s = T ge + T eg , T ge /T eg = exp(-/3/iwo) and 
ft is the inverse temperature of the environment. We use 



the trial solutions p gg = p g °J +5p g 1 g where p y gg — T ge /Ts 
is the stationary solution and Sp^J is the correction due 
to the drive. Plugging p gg in the initial equation we ob- 



(i) 



JO) 



tain 5p g X g 



PffsV^E- When integrated the adiabatic 



contribution hp gg uoq gives no contribution (see SM) and, 
with an integration by parts, the non-vanishing contri- 
bution in the limit fiEc 3> 1 can be written as 



FIG. 3: (Color online) (a) Numerically calculated (solid lines) 
dissipated average heat during the sweep in an open CPB as 
a function of the sweep time when the system is initially in 
the ground state. The system-environment coupling constant 
is chosen to be C s /Ce = 0.05 here, Ec/Ub = 1 K, and the 
environment resistance is R = 1 • 10 4 Q. The different curves 
from top to bottom correspond to e = 0.05, 0.0375, and 0.025. 
The dashed line is the analytic approximation of Eq. ( |10[ | . (b) 
Dissipated average heat at EcT/fo = 150, 100 and 50 from top 
to bottom as a function of e 2 . The other parameters and the 
line conventions are as in (a). 
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(ii) 



Here the second line arises from the detailed balance of 
the transition rates. Thus, we recover the expected prop- 
erties of the released heat: (i) it depends on the full evo- 
lution which, in this limit, is represented by the driving 
parameter q, (ii) it scales as 1/7" and (Hi) it is positive 



by a simple analytical approximation 



Q/(W) 



R C, 
Rq Cs 



E C T 



(10) 



which is derived in the SM with the assumption that, 
again, LZ-transition occurs exactly at t = T/2 with prob- 
ability Plz- Here, Rq = h/e 2 . The prediction of Eq. 
( 10 1 is shown by dashed lines in Fig. [3j in close agree- 
ment with the full numerical solution. Energy relaxation 
occurs uniformly over the positive values of q leading 
to proportionality of Q/(W) on T- As a by-product, 
Fig. [3] (a) justifies the semi-quantitative analysis and 
the proposed measurement protocol above, since most of 
the work remains stored in the system during fast ramps 
(small T). 

In the limit of fast relaxation, i.e., r -C T, it is possi- 
ble to obtain an analytical estimate of the heat released 
during the evolution. Since the dephasing time is usu- 
ally smaller than the relaxation time, we can assume 
that the coherences between the ground and excited state 
are quickly lost and they do not influence the dynamics. 



In summary, we have analyzed work done by a driving 
field on a quantum system. The obtained expression of 
average work has a physical interpretation allowing one 
to assign separate contributions to the change in the in- 
ternal energy and the heat dissipated to the environment 
in the spirit of the first law of thermodynamics. We ap- 
plied our results to a two-level system obeying in the first 
case unitary evolution and then in the presence of weak 
dissipation. For an open system, we presented a detailed 
analysis of the released heat two regimes where the re- 
laxation time was either of the order or smaller than the 
driving time. In the latter case, our approach allows an 
analytical calculation of the released heat which has and 
immediate physical interpretation. 
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Dynamical evolution in the instantaneous Landau-Zener transition model 

To make the Landau-Zener problem analytically manageable, we assume that the transition occurs instantaneously 
at time t = T/2 [TJ [2J. In this case, the unitary evolution operator can be written as U(t) = U\(t) for t < T/2, 
and as U(t) = U2(t)NLzUi(T /2) after the LZ transition at t > T/2. If the system is initially in the state \ipo) — 
a\g(0)) + /3exp(z7)|e(0)), we obtain 

\m) = e-^a\g(t)) + f3e^+^\e(t)) (12) 

for t < T /2, and 

\m) = (ay/1 - P LZ e-^+^+^ - -®)\g{t)) 

+ (av/iW^? + /Vl - PLze^f+t^Mt)) (13) 

for T/2<t< T. 

If the Cooper pair sluice is in the charging regime, i.e., Ec/Ej <C 1, the rj function is well approximated by the 
signum function up to a correction of order e 2 : rj(t) = — 1 if t < T/2 and r/(t) = +1 if t > T/2. This implies that 
the ground and the excited states can be approximated by the charge states as {\g{t)), |e(t)}} = {|0), |1)} for t < T/2 
and {\g(t}), |e(£))} = {|1), |0)} for t > T/2. This allows us to calculate directly the quantities of interest as 

f 1 - a 2 for < t < % 

<*»(*)> - <^ , , r . 2 (14) 

2a 2 )P LZ + a(a - 2^1^^008(7 + <p + ^ - P lz )Plz) for | < t < T 

Since (n(t)) does not depends on time the integration is trivial and we obtain immediately the results of Eq. (8) in 
the main text. 
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For the two-point correlator, it is convenient use the equivalent definition jj' dt 2 dt\ (n H (t2)n H (ti)) — 
2 fj~dtz J ' 2 dti$le((h H {t-^n 11 (ii))) which is obtained simply using the symmetric property of the integral and the 
fact that {{h H \t 2 )h H (ti))) f = {h H {t 2 )h H (h)) . Then, we write 



^eC(^(t 2 )ri^C*i)>) = 3FJe((Vo|^ t (i 2 )|l>(l|^(t2)C> + (ii)|l>(l|^(ti)|V'o» 

which can be calculated in a similar way. 

It turns out that, $le((h H (<2)n (*i))) does not depend on time and 



(15) 



T 



dh 



T / 

IT 
H 

4 
8 



1 



(1 - a 2 ) P LZ - a^/(a 2 - 1) (P LZ - 1) P L z cos (7 + 2ft" + <p 
(1 - 2a 2 ) P LZ +a(a- 2^(1 - a 2 ) (1 - P LZ ) P LZ cos (7 + 2^ 



(16) 

where the three contributions come from the integration intervals 0<ti<t 2 <^-,0<ti<^j<t2<T and 
< ? < ti < < T, respectively. From Eqs. (6) and (8) in the main text, the second moment is obtained as 



(W 2 ) = E 2 C P LZ . 



(17) 



Behavior of the heat in the weak coupling regime 

The full master equation that we solve numerically reads [3] 

f)gg = ~2VgJke(p ge ) ~ (T ge + T eg )p gg + T eg + T Q ^e(p ge ) 

Pge = V ge (2p gg - 1) + t^oPge ~ K T ge + ^eg)^(Pge) ~ ^ V Pge + (f + + ^-)Pgg ~ f-| 



Here 



+7^ [(T eg - T ge ) - 2(r+ + T_)p gg + 2r+ + T v (2p gg - 1) + 2(f o - f + - f_)He(p fl0 )] . (18) 
w 



9 (19) 



2g 2 + e 



is the driving term, and the various rates related to the interaction with the environment read T ge — ■w-S'v - (— Wo)j 
T eg = ^S Vg (+uJo), T v = 2^-S Vg (0), f± = ^5y g (± Wo ), f = 2^S Vg (0), T± = ^Sy g (±w ), and T - 

2 , 

(0). The couplings are defined as mi = —rjeC g /Cs and 771,2 = yl — ti^eCg/C^., and S , y B (u;) is the noise 
spectrum. In the numerical solution we have assumed temperature to be zero. With r\ = qj ' \J q 2 + e 2 , S(+ujq) = 
2RhtjQ, and 



Hloq = 2E c ^q 2 + e 2 , (20) 



we obtain the relaxation rate as 



1 , (;: fl ) 2 , 1^1.) 

where i?Q = h/e 2 . 

We consider the case in which the system is initially in the ground state and p gg (0) = 1. For the analytic 
approximation, we assume that at t — T/2 when the system passes the degeneracy, the population of the excited 
state p ee — 1 — p gg jumps from to Plz- (The ramp starts at t — and ends at t = T.) After t = T/2, the excited 
state population relaxes approximately as 

Pee I" eg Peei (^^) 

which yields 

p ee (t) = p ee (0)e-^ r - (T)tir ~p ee (0)[l- / T eg (r)dr], (23) 

Jt/2 



7 



where in the last step we have assumed that the relaxation is weak on the time scale of the sweep, T eg T <C 1. The 



dissipated heat in the ramp can be approximated by Q ~ J7-/2 Pgg( T )^ L °Q( T ) < ^ T - Then we have by inserting (20), (21 ) 

R , C 



and (23 1 into the expression of Q: 



a ,2 2 E cT 



and by dividing by the average work (W) ~ PlzEc, we have 



Q/(W) 



R , C„ 



E 



D V f~l I \ 771 J 



E C T 



(24) 



(25) 



which is Eq. (10) of the main text. 



Semiclassical approximation of the heat in the fast relaxation regime 

In the fast relaxation regime, r -C T, when the system is driven adiabatically we must include in the master 
equation the effect of the drive. It can be consider as a perturbation to the adiabatic dynamics [3 . This perturbation 



is given by the terms v ge /uiQ [where v ge is given in Eq. ( 19 ) ] which we assume to be small during the evolution. As 



a rough estimate of the adiabatic parameter we can use x — fr/(2T~Ej) since 2Ej is the smallest energy gap at the 
Landau-Zener crossing. 

At the zero-th order in x the master equation reads 

/»# = -rM+r8> (26) 

with pg°J = 0. The quasi-stationary (adiabatic) solution is pf g = T^J /I^ = 1/(1 + e~P hu} °) if we assume that the 
detailed balance T ge /T eg = exp (—fthuio) is instantaneously satisfied. The heat contribution reads explicitly 

Then we need to calculate the correction to the first order in X- This can be done starting from the master equation 
written in the superadiabatic basis {Iff*- 1 - 1 ), | e*- 1 - 1 } } [3] which allows us to include directly the high order correction in \ 
in this basis. Under the hypothesis of fast dephasing we have 

^-rM+r« (28) 

Formally we can expand the functions in the above equation as 

rgg rgg ' u Hgg 



rg> = r£> + ar«. (29) 



Using the explicit expression for the p g °J in Eq. ( 28 ) and taking the terms of the same order in \ we obtain 



^ = -4%«-^^ + ^ (30) 

which gives 



•(0) m(l) r (0) cp(l) 

°Pgg ~ r (o) j r (o)p + p(0) \ ox J 

This equation is at the first order in the adiabatic parameter x and the remaining task is to calculate the perturbation 
in the rates ST) 1 '. We have [I] 

r^J^plsU") (32, 
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where A is the noise operator acting on the system, S(uj^) is the spectral density function and lu^' is the energy 
renormalized for the drive contribution. By direct calculation and at the first order in the adiabatic parameter x, we 



,(!) 



have that 

noise is the charge noise which is described by A = |1)(1| — |0)(0|. Expanding Eq. (32) in x, we can check directly 



wo, \g w ) = \g) -iv ge /u} \e) and \gW) = 



iv ge /u>o\g}. For the Cooper pair box the dominant 



that there is no correction to the transition rates, i.e., 5T 



(i) 



S J1) = _PS9 
u rgg 



0, and, finally we obtain 



(0) 



(0) 



r 



as in the main text. With this result the released heat explicitly reads Q ■ 
parts, exploiting the detailed balance condition and the symmetry wq(T) 



(33) 



P h 2 LU (t)uJo(t) 



2r s cosh 2 (^) 



-io T ^oi[r^;§(^)]- Integrating by 
wo(0) = Ec, the first contribution reads 

(34) 



which vanishes exponentially in the limit /3Eq 3> 1. Thus, the only remaining contribution is the one in Eq. (11) of 
the main text. 
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